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ABSTRACT 


This  study  is  devoted  primarily  to  solid-solid  impact  at  velocities  in 
the  meteoric  range.  A  continuous  Eulerian  hydrodynamic  code  is  used  to 
numerically  integrate  the  flow  equations  for  several  series  of  impact  prob¬ 
lems.  The  impact  velocity,  the  projectile  shape,  the  density  of  the  pro¬ 
jectile,  and  the  equation  of  state  of  the  interacting  materials  are  varied. 
Impact  crater  dimensions  are  found  to  increase  as  the  0.  58  power  of  the 
impact  velocity.  Projectiles  of  the  same  mass  but  different  densities  have 
the  same  overall  impact  effects  on  the  target,  for  velocities  above  about 
10^  cm/  sec. 

One-dimensional  and  axisymmetric  impacts  are  also  computed  for 
materials  with  an  ideal-gas  equation  of  state.  These  calculations  provide 
some  important  comparisons  with  analytical  work  on  the  problem  of  hyper¬ 
velocity  impact. 

Many  of  the  concepts  used  in  solid- solid  impact  are  already  present 
in  the  one -dimensional  problem.  The  theory  for  this  well-understood 
impact  is  discussed  in  considerable  detail. 
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1.  INTRODUCTION  AND  SUMMARY 


The  high-velocity  impact  interaction  is  separable  into  an  early  phase 
which  is  governed  by  the  equations  of  fluid  dynamics  and  a  later  phase 
during  which  material  strength,  not  contained  in  the  fluid  dynamics,  is  no 
longer  a  negligible  effect. 

The  present  report  is  devoted  primarily  to  the  early  fluid- dynamic , 
or  hydrodynamic,  phase  of  the  impact  interaction.  Special  emphasis  is 
placed  on  the  property  of  "late-stage  equivalence"  within  this  hydrodynamic 
phase,  whereby  it  is  found  that  projectiles  of  differing  mass  and  impact 
velocity  or  of  various  shapes  or  of  unequal  densities  can  give  rise  to 
target-plate  flows  which  are  very  nearly  identical  after  early  transient 
differences  have  died  out. 

The  property  of  late-stage  equivalence,  when  it  occurs  within  the 
hydrodynamic  phase  of  the  interaction,  is  useful  because  it  enables  one  to 
circumvent  the  problem  of  strength-dependent  deformation  and  to  establish 
the  dependence  upon  the  initial  conditions  of  overall  impact  effects,  such 
as  total  impulse  to  the  target  or  target  crater  size.  For  example,  it  is 
possible  to  show  that  an  impact  in  the  meteoric  velocity  range  has  the 
same  overall  effects  as  a  suitable  impact  at  velocities  already  attained  in 
controlled  experiments.  Any  alternative  method  of  learning  the  overall 
effects  of  impacts  at  meteoric  velocities  would  involve  either  treating  the 
strength-dependent  deformation  problem  theoretically  or  solving  the  formi¬ 
dable  experimental  problem  of  accele  '.ting  solid  projectiles  to  these 
velocities . 

The  governing  equations  and  the  solid  equations  of  state  used  in  this 
study  are  given  in  Sections  2  and  3.  In  Section  4  the  continuous 
Eulerian  hydrodynamic  code  is  described  and  results  for  a  typical 
impact  problem  are  compared  with  those  from  the  same  problem 
done  with  the  particle  -  in- cell  code.  The  continuous  Eulerian  version  gives 
a  substantially  smoother  solution  which  is  in  good  general  agreement  with 
integrated  or  smoothed  results  from  the  particle-in-cell  code.  Furthermore, 
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the  Eulerian  code  requires  considerably  less  computing  time  and  has  been 
used  in  almost  all  of  the  impact  work  discussed  in  this  report. 

One-dimensional  and  axisymmetric  impacts  of  materials  with  ideal- 
gas  equation  of  state,  PV  =  (y  -  1)E,  are  discussed  in  Sections  5  and  6, 
respectively.  The  assumption  that  impacts  are  late-stage  equivalent  on 

'Vy 

the  basis  of  equal  MqVq  ,  rather  than  dependent  at  late  times  on  the  pro¬ 
jectile  mass,  Mp,  and  impact  velocity,  v^,  individually,  is  shown  to  lead 
to  the  result  that  the  flow  is  of  the  familiar  self-similar  type.  For  slab 
impact  and  for  the  special  case  y  =  1.4,  this  self-similiar  flow  has  been 
treated  analytically  in  recent  years,  so  that  a  comparison  can  be  made 
between  the  analytical  solution  and  that  obtained  numerically  without  the 
assumption  of  late-stage  equivalence.  The  numerical  solution  and  the 

analytical  similarity  solution  are  in  excellent  agreement  at  late  times, 
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and  impacts  are  found  to  be  equivalent  on  the  basis  of  equal  M^v^'  ,  an 

exponent  which  is  midway  between  those  for  equivalence  on  the  basis  of 
equal  kinetic  energy  and  equal  momentum.  Additional  considerations 
indicate  that  this  exponent  varies  from  1.0  to  1.790  as  y  is  varied  from 
1 . 0  to  infinity . 

For  axisymmetric  ideal-gas  impact,  the  similarity  solution  apparently 
cannot  be  obtained  analytically  without  making  additional  simplifying  approx¬ 
imations.  But  it  is  possible,  from  the  form  of  the  solution  indicated  by 
the  assumption  of  late-stage  equivalence,  to  determine  the  relation  which 
must  exist  between  the  shock  attenuation  rate  and  a  and  between  the  rate 
of  increase  of  selected  integrated  momenta  and  a.  These  results  can  be 
compared  with  the  computed  solution,  determined  without  assuming  late- 
stage  equivalence.  Good  agreement  is  found  at  late  times  with  the  self¬ 
similar  property  of  the  similarity  solution,  and  the  value  of  a  determined 
from  the  shock  attenuation  and  momentum  data  is  a  =  0.  59.  This  a  is  to 
be  compared  with  the  value  1/3  tor  momentum  equivalence  and  2/3  for 
energy  equivalence.  These  results  are  for  y  =  1.5  and  should  therefore 
apply  to  solids  at  extremely  high  impact  velocities  because  of  equation- 
of-state  similarities. 

In  the  case  of  solid-solid  impact,  discussed  in  Section  7,  the  assump- 

3/y 

tion  of  late-stage  equivalence  on  the  basis  of  equal  M^v^  does  not  lead  to 
the  requirement  that  the  solution  be  expressible  as  a  similarity  solution. 
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The  assumption  itself,  however,  is  found  to  be  valid  at  high  velocities  by 
comparing  computed  solutions,  and  this  result  can  be  used  to  establish 
overall  impact  effects.  The  results  of  computations  for  some  forty  impact 
interactions  are  summarized  in  Section  7.  Particular  attention  is  devoted 
to  establishing  late -stage  equivalence  and  the  limits  of  its  validity,  not 
only  in  regard  to  impact  velocity  variation,  but  also  with  respect  to  changes 
in  projectile  shape  and  equation  of  state  of  the  impacting  solids. 

In  Section  8  the  hydrodynamic  equivalence  results  from  Section  7  are 
used  to  determine  relations  for  such  overall  effects  as  target  crater  size, 
total  impulse  to  the  target,  and  the  effects  of  target  hardness  on  crater 
size.  Some  concluding  remarks  are  given  in  Section  9. 

The  present  effort  is  an  extension  of  the  work  reported  in  Ref.  1. 

This  initial  study,  which  was  based  primarily  on  the  application  of  a 
particle-in-cell  code  to  iron-iron  impact,  is  in  good  general  agreement 
with  results  given  here.  The  work  reported  in  the  present  document  has 
been  devoted  to  two  main  objectives --namely,  the  development  of  a  more 
accurate  and  economical  (Eulerian)  hydrodynamic  code  for  the  treatment 
of  solid- solid  impact,  and  the  application  of  the  code  to  a  systematic  study 
of  the  more  important  parameters  defining  an  impact  interaction.  The 
more  significant  particle -in-cell  impact  problems  have  been  redone  in 
order  to  check  the  earlier  results  and  the  investigations  are  extended  to 
other  materials,  to  varied  projectile  shapes,  and  to  unlike-material 
impact. 

It  is  a  pleasure  to  acknowledge  the  aid  of  Dr.  Floyd  E.  Allison  of  the 
Ballistic  Research  Laboratories,  who  is  the  technical  monitor  for  this  study. 
We  are  grateful  to  Dr.  Allison  for  a  number  of  very  helpful  discussions 
of  pertinent  theoretical  and  experimental  work  on  the  subject  of  hyper¬ 
velocity  impact. 
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intensive  quantities  such  as  pressure,  density,  and  velocity  remain  unchanged. 
This  scaling  property  is  used  in  subsequent  sections  to  change  the  relative 
size  of  two  impact  interactions  so  that  flow  configurations  are  in  agreement 
at  late  times. 

The  principal  approximations  used  to  derive  the  governing  equations  are 
that  viscosity,  heat-conduction,  radiation,  and  strength  effects  within  the 
continuous  flow  are  neglected.  Neglecting' viscousand  heat -conduction  terms 
can  be  justified  for  flows  of  sufficiently  large  dimensions.  Quantitative  con¬ 
siderations^  suggest  that  the  effects  are  negligible  for  flows  of  dimensions 
-4  -3 

greater  than  10  or  10  cm,  although  there  exists  considerable  uncertainty 

in  this  lower  limit  because  the  appropriate  viscous  coefficient  for  metals  is 

very  uncertain.  The  transport  by  radiation  of  a  significant  fraction  of  the 

energy  within  the  interaction  is  possible  only  at  the  very  high  temperatures 

7  7 

associated  with  impacts  around  3  x  10  to  5  x  10  cm/sec  or  greater.  The 
exclusion  of  material  strength  from  the  description  is  possible  so  long  as 
the  stresses  are  very  large  compared  to  the  yield  strength  of  the  solid. 

This  approximation  is  valid  for  early  stages  of  the  interaction,  but  fails 
at  late  times.  The  present  treatments  of  the  flow  are  limited  to  this  early 
hydrodynamic  phase. 

The  above  approximations  and  the  simple  scaling  law  are  closely  re¬ 
lated.  Viscous,  heat- conduction,  and  radiation  effects  would  introduce 
higher-order  derivatives  in  the  continuous -flow  equations,  causing  depar¬ 
tures  from  the  scaling  law.  Although  strength  effects  are  not  important 
in  the  hydrodynamic  phase,  it  is  pertinent  to  note  (in  connection  with  later 
estimates  of  overall  impact  effects)  that  ordinary  stress-strain  effects  do 
scale,  whereas  time  dependence  within  the  stress-strain  relation  causes  a 
breakdown  of  simple  scaling.  Thus  the  experimental  determination  that 
impact  effects  scale,  e.g.,  that  crater  size  is  directly  proportional  to 
projectile  size  or  that  target  impulse  varies  directly  with  projectile  mass, 
suggests  that  the  neglected  effects  do  not  play  an  important  role  in  the 
interaction.  Observed  departures  from  such  scaling  are  very  important, 
on  the  other  hand,  as  evidence  that  new  mechanisms,  such  as  viscous 
forces  or  rate-sensitive  strength  effects,  must  be  included  in  the  equations. 
The  main  portion  of  experimental  evidence  to  date  appears  to  indicate  tha* 
the  simple  linear  scaling  is  at  least  a  good  approximation. 


3.  EQUATION  OF  STATE 


During  the  impact  inte  raction,  each  element  of  the  material  is  processed 
by  a  strong  shock  wave  and  is  then  free  to  expand  isentropically.  It  is 
therefore  desirable  to  have  an  equation-of- state  representation  which  ac¬ 
curately  reproduces  both  the  Hugoniot  curve  and  expansion  isentropes  of 
the  material.  Such  equation-of- state  fits  for  ten  common  metals  and  a 
typical  plastic  (polyethylene)  have  been  made  and  described  in  considerable 
detail  in  Ref.  2.  The  input  data  for  the  fits  consist  of  available  experimen¬ 
tal  shock-wave  results,  usually  extending  to  a  few  megabars,  and  of  theo¬ 
retical  equation-of- state  results  by  Cowan  from  the  Thomas -Fermi-Dirac 
model  of  the  atom,,  for  pressures  of  several  tens  of  megabars  and  greater. 
For  condensed  states,  the  equation  of  state  is 


P  = 


a  + 


V 


+  1 


^  +  Au  +  B  yiz 


and  for  expanded  states,  the  pressure  is  given  by 

P  =  aEo  + 


au.-b'v/vo-1> 


+  1 


Le0ti‘ 


e"0!(  V/ Vq-  1) 


where  P  is  pressure,  E  is  specific  internal  energy,  p  is  density  (p  =  1/V), 
and  T|  =  p  /pg,  (j  =  71  -  1 .  The  remaining  parameters  are  chosen  to  fit  the 
equation-of-state  data  for  the  individual  materials.  Values  of  these  param¬ 
eters  for  the  eleven  materials  are  given  in  Table  1. 

For  the  calculations  in  this  study,  the  condensed  form  of  the  equation 

of  state  has  been  used  for  states  V/V„  <  1  and  for  all  states  E  <  E  and  the 

U  s 

expanded  form  of  the  equation  of  state  is  usedforV/Vg>  landE>Eg..  This 
definition  of  regions  for  the  condensed  and  expanded  formulations,  is 
slightly  different  from  that  given  previously,  the  change  having  been  made 
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to  avoid  an  artificial  vaporization  of  highly  expanded,  but  cold,  debris. 
Finally,  negative  pressures  which  may  arise  in  the  condensed  formulation 
have  been  set  to  zero  in  the  calculations,  except  in  two  problems  where 
the  effects  of  this  simplification  were  explored. 


Table  1 

EQUATION -OF-STATE  CONSTANTS 


Material 

a 

■ 

H 

E0  3 

(Mb-cm  /g) 

A 

(Mb) 

B 

(Mb) 

E 

s 

(Mb-cmVg) 

o t 

3 

ch2 

EQ 

i 

2.0 

0.92 

VI 

0.075 

0.02 

0.  024 

10 

5 

Pb 

09 

2.  37 

1 1 .  34 

\-ma  1 

0.4664 

0.15 

0. 0026 

13 

15 

w 

0.  5 

1.04 

19.  17 

1 

3.08 

2.50 

0.014 

10 

10 

Cu 

0.5 

1.5 

8. 90 

BBS  i 

1.39  i 

1.  10 

0.  02 

5 

5 

Fe 

0.  5 

1.5 

7.86 

0.095 

1.279 

1.05 

0.025 

5 

5 

A1 

■ 

1.63 

2.  70 

0.050 

0.  752 

0.65 

0.03 

5 

5 

Be 

Hi 

0.62 

1.85 

0.  175 

1.1734 

0.55 

0.  10 

5 

5 

Ti 

09 

0.60 

4.  50 

0.070 

1.03 

0.50 

0.035 

5 

5 

Ni 

0.  5 

8.86 

0.090 

1.912 

1. 50 

0.03 

5 

5 

Mo 

□ 

10.  20 

0.045 

2.713 

1.65 

0.03 

5 

5 

Th 

Bi 

11.68 

0.025 

0.531 

0.50 

0.02 

9 

0.  88 

4.  HYDRODYNAMIC  CODES 


An  excellent  general  description  of  the  particle -in-cell  method  is  to 

(3  4) 

be  found  in  the  original  Los  Alamos  Scientific  Laboratory  literature.  '  ' 

Two  similar  particle -in-cell  codes,  SHELL  and  SPEAR,  are  operative  at 

General  Atomic,  and  the  continuous  Eulerian  code  uses,  insofar  as  possible, 

( 5) 

the  general  scheme  found  in  SHELL. 

A  central  notion  in  all  of  these  codes  is  the  division  of  space  into  cells 
which  remain  fixed  and  through  which  the  fluid  moves.  Each  cell,  at  time 
t,  has  associated  with  it  values  of  the  fluid  density,  total  fluid  energy,  and 
the  components  of  velocity  of  the  fluid.  The  process  of  advancing  the  de¬ 
pendent  variables  to  time  t  +  At  is  done  in  two  distinct  phases .  In  PhaseT,  the 
field  terms  in  the  equations,  all  involving  pressure,  are  used  to  partially 
update  the  dependent  variables,  with  the  remaining  transport  terms  being 
neglected.  In  Phase  II,  the  transport  terms  are  accounted  for  in  a  process 
which  involves  moving  material  across  cell  boundaries.  Finally,  new  cell 
pressures  are  computed  from  the  updated  densities  and  energies  and  the 
time  step  for  the  next  cycle  is  determined. 

The  continuous  Eulerian  and  the  particle-in-cell  codes  differ  in  the 
procedure  used  for  the  Phase  II  material  transport.  In  the  particle-in¬ 
cell  scheme,  each  particle  represents  a  finite  mass  so  that  very  small 
variations  in  cell  density  are  not  accurately  reproduced.  The  solutions 
obtained  with  the  particle -in-cell  code  exhibit  unreal  cell-to-cell  fluctua¬ 
tions,  but  integrated  or  suitably  averaged  quantities  appear  to  be  reliable. 
The  continuous  Eulerian  code,  on  the  other  hand,  permits  small  variations 
in  cell  density  and  rarefaction  to  very  small  densities.  A  disadvantage  to 
the  Eulerian  scheme  is  the  loss  of  the  Lagrangian  aspects  of  the  particle 
scheme.  Elements  of  mass  are  no  longer  individually  identified  from 
cycle  to  cycle,  and  the  location  of  material  interfaces  and  material  free 
surfaces  also  offers  special  problems.  In  the  present  effort,  special 
restrictions  on  the  transport  have  been  used  to  prevent  the  artificial 
smearing  out  of  solid  free  surfaces,  and  applications  of  the  Eulerian  code 
have  so  far  been  limited  to  like -material  impacts. 
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A  number  of  problems  have  been  computed  in  the  process  of  developing 
the  continuous  Eulerian  code,  and  the  results  have  been  compared  with 
available  analytical  solution  and  with  computed  results  from  the  SHELL, 
and  SPUTTER  (one -dimensional  Lagrangian)  hydrodynamic  codes.  Flows 
treated  include  the  shock-tube  problem,  the  one -dimensional  rarefaction 
to  zero  pressure,  the  spherical  blast  problem,  and  a  number  of  plane  and 
axisymmetric  impact  interactions.  In  general,  these  comparisons  have 
been  quite  encouraging,  with  good  agreement  being  found  between  the  con¬ 
tinuous  Eulerian  results  and  the  analytical  or  SPUTTER  solutions.  The 
particle-in-cell  and  continuous  Eulerian  results  are  also  in  good  general 
agreement,  with  the  former  solution  being  relatively  free  from  undesirable 
oscillations.  The  most  extensive  comparison  of  the  particle-in-cell  and 
continuous  Eulerian  codes  has  been  for  the  axisymmetric  impact  problem 
described  in  the  following  paragraphs. 

An  iron  projectile,  a  right  circular  cylinder  3  cm  long  by  3  cm  diameter, 
impinges  on  a  semi-infinite  iron  target  at  a  velocity  of  4  x  10^cm/sec.  The 
problem  was  run  using  both  the  continuous  Eulerian  and  the  particle-in-cell 
codes.  The  grid  size  and  starting  conditions  we^e  identical  for  both  prob¬ 
lems,  with  72  cells  in  the  pellet  and  2976  ceils  in  the  target.  Cell  dimen¬ 
sions  were  0.25  cm  and  the  particle  code  used  16  particles/cell. 

Representative  comparisons  of  results  from  the  continuous  Eulerian 
and  particle  codes  are  given  in  Figs.  1,  2,  and  3.  Figure  1  is  a  plot  of 
momentum  versus  time  data  for  the  two  problems.  Here  the  axial  momen¬ 
tum  £  M-v^,  is  the  momentum  in  the  original  pellet  flight  direction  summed 
over  all  material  with  a  positive  component  of  velocity  in  that  direction. 
Similarly,  the  radial  momentum  £  IvLu^  is  an  integrated  momentum  for  the 
radial  motion  and  is  summed  over  all  material  above  the  initial  target- 
projectile  interface  and  for  which  the  radial  velocity  is  positive.  The 
agreement  between  momentum  results  from  the  two  codes  is  seen  to  be 
remarkably  good.* 

Figures  2  and  3  are  plots  of  pressure,  compression,  and  velocity 
profiles  from  the  two  codes  at  successive  times  and  along  a  column  of 


5{C 

These  momenta  have  proved  useful  as  accurate  measures  of  the 
total  target  disturbance  and  are  used  elsewhere  in  this  report. 


PARTICLE  IN  CELL 


-0-0  0-  CONTINUOUS  EULERIAN 


TIME  (/iSECt 


Comparison  of  radial  and  axial  momenta  from  the  particle -in-cell 
and  continuous  Eulerian  codes 


r 


Fig.  ^.-Comparison  of  axial  velocity,  compression,  and  pressure  profiles  from  the  particle-in-cell 

and  continuous  Eulerian  codes 


(SMV0W)  3W0SS38d 
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and  continuous  Eulerian  codes 


13 


cells  near  the  axis  of  symmetry.  The  results  given  here  are  typical  of 
those  which  have  been  obtained  in  a  number  of  such  comparisons:  The 
continuous  Eulerian  solution  is  relatively  smooth  and  in  good  general 
agreement  with  the  particle-in-cell  results,  the  cell-to-cell  fluctuations 
in  the  latter  being  attributed  to  the  approximate  mass  resolution.  The 
fluctuations  are  more  pronounced  at  the  later  time  owing  to  the  necessity 
of  defining  smaller  compressions. 

The  particle-in-cell  code  in  the  above  comparison  required  longer 
machine  time  by  a  factor  of  15,  the  additional  time  being  consumed  in 
processing  data  for  the  individual  particles. 

A  more  detailed  description  of  the  continuous  Eulerian  code  is  planned 
by  W.  E.  Johnson. 


5.  ONE -DIMENSIONAL  IMPACT  WITH  IDEAL-GAS 
EQUATION  OF  STATE* 


The  flow  resulting  from  the  impact  of  a  slab  on  a  semi-infinite  target 
contains  many  of  the  essential  elements  of  three-dimensional  impact  theory 
and  yet  is  amenable  to  a  relatively  simple  analytical  treatment  provided 
the  target  and  slab  are  both  taken  to  have  an  ideal-gas  equation  of  state 
and  to  be  devoid  of  internal  energy  prior  to  impact.  In  addition,  it  is  pos¬ 
sible  to  obtain  very  accurate  solutions  to  the  flow  equations  by  numerical 
integration  using  SPUTTER  or  other  hydrodynamic  codes. 

The  first  portion  of  this  section  contains  the  formulation  of  the  equa¬ 
tions  for  plane  similarity  flow  based  on  the  notion  of  late-stage  equivalence. 

OL 

The  fact  that  only  one  projectile  parameter,  LqVq,  is  allowed  to  appear  in 
the  problem  is  sufficient  to  reduce  the  general  flow  equations  from  a  set 
of  partial  differential  equations  to  a  first-order  ordinary  differential  equa¬ 
tion.  The  value  of  a  is  determined  by  solving  this  equation  with  appropriate 
boundary  conditions. 

The  second  portion  of  this  section  makes  use  of  the  previous  mathe¬ 
matical  results  to  examine  the  mechanics  of  the  flow.  It  is  shown  that  the 
flow  has  infinite  energy  and  zero  momentum,  but  that  the  product  of  energy 
times  momentum  to  a  certain  power,  El  ,  is  finite.  This  product  is  a 
function  of  L^v®.  In  one  case  where  the  equations  can  be  solved  exactly, 
d  is  found  to  be  3/2.  In  this  case,  it  is  shown  that  it  is  correct  to  equate 
the  El  product  for  the  similarity  flow  to  the  value  of  E^I^  of  the  slab.  In 
general,  a  depends  on  the  equation  of  state  chosen.  For  a  perfect  gas, 

0/  is  obtained  as  a  function  of  y  by  interpreting  the  work  of  previous  authors 
in  terms  of  the  concept  of  late-stage  equivalence. 


>!< 

This  section  is  essentially  self-contained.  It  can  be  read  separately 
from  the  rest  of  the  report  and  can  be  skipped  in  a  first  reading  of  the 
report. 
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5.1.  MATHEMATICAL  ANALYSIS  OF  SIMILARITY  FLOW 

The  analysis  of  the  flow  is  most  conveniently  discussed  by  formulating 
the  Lagrangian  equations  of  motion.  Let  h  denote  the  number  of  grams  per 
square  centimeter  separating  a  particle  from  the  free  surface  at  h  =  0, 
u(h,  t)  be  the  velocity  of  the  particle  specified  by  h,  p  is  the  pressure,  V 
is  the  specific  volume,  and  its  initial  value.  It  follows  that  the  initial 
distance  of  a  particle  from  the  free  surface  is  hV-.  The  equations  of 
motion  can  then  be  written:' 


v  -  “h 

(conservation  of  mass)  , 

(1) 

u  =  -ph 

(conservation  of  momentum)  , 

(2) 

p  =  f(h)V"y 

(isentropic  expansion)  . 

(3) 

The  boundary  conditions  for  the  flow  at  the  shock  front  are  the  Rankine- 
Hugoniot  relations.  These  are  reproduced  below  for  a  target  that- is  a  per¬ 
fect  gas  with  no  initial  internal  energy: 


The  subscript  s  denotes  the  value  of  any  quantity  at  the  shock  front  and  vg 
is  the  shock  velocity.  In  addition,  the  pressure  must  vanish  on  the  free 
surface. 

The  independent  variables  of  the  problem  are  h  and  t,  but  the  principle 
of  late -stage  equivalence  leads  to  the  conclusion  that  h  can  only  appear 
with  Vq,  Lq,  and  v^  in  the  nondimensional  combination 


L  a  t 

°V°  (hv0)a+1 


It  is  natural  to  choose  as  a  similarity  variable,  r),  the  ratio  of  h  to  its  value 
at  the  shock  front,  hg(t).  This  approach  can  be  reconciled  with  late-stage 
equivalence  by  determining  hg(t)  from  the  equation 
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h 


1 

w 


vP* 


n-l/(a+l) 


<hV 


Of+1 


(7) 


where  0  is  an  undetermined  constant.  The  position  of  the  shock  is  then 
given  explicitly  by 


hs(t)  =  0(Lovot  J  /v o  • 


(8) 


If  a  similarity  solution  exists  and  late-stage  equivalence  holds,  then 
the  pressure,  velocity,  and  specific  volume  must  have  the  form 


a  2/a 

_  e  p(w) 

p  -  81  v  l+2  Int  ’ 

v  0  ns 

a  1/a 

(L0V0)  u(n) 

-  *2  v‘/«  ^ 


(9) 


(10) 


V  =  03VoV(t?) 


(ID 


where  0j,  and  0^  are  constants  and  the  barred  variables  are  non- 

dimensional  functions  defined  to  be  unity  at  the  shock  front,  i.e.  , 

F(i)  =  TKD  =  V(D  -  l  ..  (12) 

In  order  that  the  Rankine -Hugoniot  conditions  be  satisfied  at  the  shock,  it 
is  necessary  that 

\a/(cm  1) 


9 

0, 


V  ±  1  «  +  1  q  I 
2  a  “zj 


-  y  +  1  02 
2  W2  ’ 


(13) 

(14) 


0  =  V  -  1 

°3  y  4-1 


(15) 
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The  similarity  assumption  inherent  in  Eqs.  (9),  (10),  and  (11)  reduces 
the  flow  equations  to  the  following  second-order  system  of  ordinary  differ¬ 
ential  equations  for  p,  u,  and  V: 

T)tP  +  xt/a  =  V',  (16) 

u’  +l-iJvV'  =  0  ,  (17) 

vz/avvy  =  i  .  (is) 

These  equations  can  be  solved  for  the  derivatives  p' ,  tP ,  and  V': 


P' 


_  £  £ 
a  v 


1  + 


1  - 


y  u 
y  -  i  r/fr 
2y 

y  -  1 


~zZ 
V  v 


(19) 


V' 


I  2 

a  y  -  1 


y  -  1 


(20) 


u' 


Vi 

an 


1 

JE. 


i 


-b 

V  V 


(21) 


Inspection  of  these  equations  suggests  introducing  a  new  set  of  variables 
q  and  y,  whe  re 


q 


7  -2 

Zy  u 

y  - 1  v  v 


z 

y  . 


(2Z) 


y  =  r?u/p  , 


(23) 


and  q  is  proportional  to  the  square  of  the  local  Mach  number,  M.  After 
some  calculation,  one  finds  the  following  first-order  equation  for  q(y): 


=  a  -  4y  +  y  '  ^)q  . 

dy  y  (1  +  a)y2  -  Zy  +  (y  -  a)q 


(24) 
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The  choice  of  variables  given  by  Eqs.  (22)  and  (23)  was  suggested  by 
M.  Nowak^^  and  results  in  a  somewhat  simpler  equation  than  given  by 
previous  authors.  The  solution  must  satisfy  the  boundary  condition 


Furthermore,  the  pressure  must  vanish  at  the  free  surface.  In  terms  of 
the  present  variables,  this  requires  that  the  behavior  of  the  equations  must 
be  determined  when  the  magnitudes  of  q  and  y  are  large.  It  can  be  shown 
asymptotically  that 

-  (1  +a)y  .  (26) 

Combining  this  result  with  Eqs.  (19)  and  (21)  provides  the  required  be¬ 
havior  of  the  pressure  and  velocity  at  the  free  surface: 


P  ~  V 


(2/a)[(y/2)-l]+y 


tr  ~  t) 


This  result  will  be  used  below. 


Equation  (24)  has  a  saddle  point  at 


q  =  4  , 


and  it  can  be  shown  that  the  solution  curve  must  pass  through  this  singu¬ 
larity.  If  it  does,  the  pressure  vanishes  at  the  free  surface  and  the  flow 
is  everywhere  well  behaved.  Thus,  the  free-surface  condition  can  be  re¬ 
placed  by  Eq.  (29).  Equations  (25)  and  (29)  together  constitute  two  bound¬ 
ary  conditions  for  the  first-order  system  represented  by  Eq.  (24).  It  is 
possible  to  satisfy  both  boundary  conditions  if  a  is  appropriately  chosen 

and  a  is  thereby  determined.  In  general,  the  problem  cannot  be  solved 

( 8) 

analytically,  but  numerical  approaches  have  been  used  by  H&fele'  and 

(9) 

by  Rae  and  Kirchner'  '  to  solve  equations  equivalent  to  Eq.  (24).  In  their 
analysis  a  similarity  approach  based  on  the  parameter  x/tn  was  used. 
Their  approach  is  consistent  with  the  concept  of  late-stage  equivalence. 

A  comparison  of  the  similarity  parameter  x/tn  with  Eq.  (7)  shows  that 


a 


n 

-  n 


(30) 


5.2.  PHYSICAL  INTERPRETATION  OF  SIMILARITY  FLOW 

It  is  possible  to  plot  a  as  a  function  of  y,  using  the  detailed  tables  of 

/  9) 

Rae  and  Kirchner.  The  result  is  given  in  Fig.  4.  It  is  found  that  for 
y  near  one,  a  is  also  near  one  and  impacts  are  nearly  equivalent  on  a 
momentum  basis,  a  increases  with  y  and  asymptotically  approaches  the 
value  1. 790,  so  that  for  no  value  of  y  are  flows  ever  exactly  equivalent 
on  an  energy  basis. 

The  total  energy  of  the  flow  to  the  right  of  h  is  given  by  the  integral 


E 


and  the  momentum  to  the  right  of  h  is  given  by 


I 


(31) 


(32) 


As  h  goes  to  zero,  these  quantities  approach  infinity  and  zero,  respectively. 
The  rate  at  which  these  limits  are  approached  is  given  by  substitution  from 
Eq.  (28): 


E  sa  h 


l-(2  fa) 

f 


(33) 


I 


hM1/a). 


(34) 


From  the  principle  of  late -stage  equivalence  it  is  expected  that  the 
flow  should  depend  only  on  energy  and  momentum  in  a  certain  combination, 
El\  where 


El 


Lovo 


<LoV 


[- 


(2+X)/( 14 X 
0  0 


( 1+X) 

>] 


(35) 


The  relation  between  X  and  a  is  thereby  determined  as 
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This  consideration  suggests  that  ei^-2) /(a-  1)  obtained  from  the  order  of 
magnitude  integrals  above  should  be  a  finite  quantity;  that  is  indeed  the  case 
as  verified  by  examining  Eqs.  (33)  and  (34). 

It  seems  plausible  that  the  product  EI^2  1)  taken  over  the  late- 

2  /  ^  \  j  /  ^  y  V 

stage  flow  field  should  be  the  same  as  the  product  1jqvq  *  (LqVq)'  '  '  '  2 

determined  from  the  energy  and  momentum  of  the  impacting  slab.  This 
cannot  be  proved  without  integrating  the  exact  flow  equations,  which,  in 
general,  cannot  be  done.  However,  it  is  possible  to  carry  through  the 
calculations  in  one  particular  case.  A  number  of  authors  have  observed 
that  the  similarity  equations  can  be  integrated  if  the  pressure  is  linear 
in  h  or,  in  terms  of  similarity  variables,  if 'p  =  17.  This  assumption  leads 
to  a  solution  consistent  with  all  the  boundary  conditions  only  if  y  =  7/5 
and  a  =  3/2,  and  if  so, 

2/3 

u  =  92(L0V0/2/V0)  ^r,)/h^/3,  (37) 


where 


u(n)  =  (3  -  T)’2/3)/2  .  (38) 

The  asymptotic  forms  of  the  energy  ar.d  momentum  integrals  are  then, 
for  tj  small, 


E  ~ 


Vo3'2| 

00 

1  v0 

h*'V/3 

s  1 


(39) 


I  ~ 


W2' 


2/3 
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-yv'3 


(40) 


Equating  El  (X  =  1  in  this  case)  to  the  energy-momentum  product  for  the 
slab  determines  9^: 


93  =  8/9 


for  y  =  7/5  . 


(41) 
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El  may  be  thought  of  as  an  integral  invariant  of  the  flow.  This  interpretation 
is  strengthened  by  considering  the  energy  E+  and  momentum  I+  of  that 

portion  of  the  material  which  is  moving  in  the  direction  of  the  original  slab. 

-1/5  1/5 

The  energy  decreases  as  t  and  the  momentum  increases  as  t  .  How¬ 

ever,  the  product  E+I+  is  constant  and  equal  to  0.  5132  of  the  product  taken 
over  the  entire  flow. 

To  verify  this  hypothesis  that  the  El  product  for  the  slab  and  the  simi¬ 
larity  flow  are  equal,  a  numerical  integration  of  the  flow  equations  was 
carried  out  using  the  SPUTTER  code.  This  has  also  provided  an  example 
of  how  rapidly  the  flow  approaches  the  similarity  solution.  The  target 
material  was  chosen  to  have  a  density  of  one;  the  velocity  of  the  projectile 
was  taken  to  be  10^  cm/ sec;  and  the  mass  per  unit  area,  hg  =  Lg/Vg,  was 
taken  to  be  1.2  g/cm^.  At  impact  the  shock  strength  was  0.  3  Mbar. 

Pressure,  velocity,  and  density  data  obtained  from  the  computer  solu¬ 
tion  are  presented  in  Figs.  5,  6,  and  7  and  are  compared  with  the  analytical 
solution  at  corresponding  times.  In  the  figures,  time  is  given  from  the 
instant  of  contact  in  the  finite -slab  solution.  The  appropriate  time  origin 
for  the  analytical  solution  is  the  instant  of  impact  by  a  slab  of  zero  thick¬ 
ness  and  is  found  by  subtracting  0. 88  jjsec  from  the  numbers  that  are 
given.  All  data  points  obtained  from  the  SPUTTER  code  solution  in  the 
neighborhood  of  the  shock  are  shown.  However,  only  representative  data 
points  are  given  where  the  oscillations  found  near  the  shock  front  have 
died  out.  No  attempt  was  made  to  optimize  the  quality  of  the  computer 
solution,  since  the  agreement  with  the  analytical  result  is  satisfactory. 
Increase  of  artificial  viscosity  in  the  SPUTTER  code  would  reduce  the 
shock-front  oscillations. 

These  results  can  be  interpreted  by  considering  a  family  of  finite- 

slab  impacts  with  different  values  of  thickness,  L,g,  and  velocity,  Vg, 

OL 

subject  to  the  constraint  that  LgVg  is  the  same  for  all  the  slabs.  At  early 
times  the  flow  profiles  are  all  quite  different,  but  at  late  times  they  all 
become  similar.  The  thin  high-velocity  slabs  result  in  flows  which  ap¬ 
proach  the  similarity  flow  soonest.  In  fact,  the  similarity  solution  can 
be  interpreted  as  an  impact  by  a  slab  of  zero  thickness  and  infinite  veloc- 
ity.  With  the  constraint  that  LgVg  is  constant,  it  also  follows  that  the 
similarity  solution  corresponds  to  impact  by  a  slab  with  infinite  energy 
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Fig.  5- -Comparison  of  pressure  in  the  finite -slab  impact  problem  with  the  similarity  solution 
for  y  =  7/5  at  three  times;  h  is  measured  from  the  free  surface 
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ab  impact  problem  with  the  similarity  solution 
measured  from  the  free  surface 
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Fig.  7- -Comparison  of  density  in  the  finite- slab  impact  problem  with  the  similarity  solution 
for  y  =  7/5  at  three  times;  h  is  measured  from  the  free  surface 
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and  zero  momentum.  Since  a  similarity  solution  can  contain  only  one  free 
dimensional  parameter,  in  this  case  L^Vq,  the  energy  and  momentum  of  the 
impacting  slab  cannot  both  be  reconstructed  from  looking  at  the  late -stage 
similarity  flow.  However,  the  product  EI^  Oi)/(d  1)  characteristic  of 
the  impact  and  reflects  the  dependence  of  the  flow  on  the  initial  momentum 
and  energy  of  the  slab. 


6.  AXISYMMETRIC  IMPACT  WITH  IDEAL-GAS 
EQUATION  OF  STATE 


i 

I 


The  impact  behavior  of  an  initially  cold  ideal-gas  material  is  a 
considerably  simpler  problem  than  that  of  solids  because  the  effect  of 
variation  in  impact  velocity  can  be  determined  from  dimensional  considera¬ 
tions  alone.  Thus,  it  is  necessary  to  solve  the  problem  for  only  one  im¬ 
pact  velocity  for  each  y  and  initial  geometry. 

For  a  given  geometry,  say  that  of  a  sphere  cr  cylindrical  pellet  on  a 
semi-infinite  plate,  the  problem  is  fully  characterized  by  given  values  of 
y,  the  initial  density,  p^,  the  size  (a  length  L^)  of  the  projectile,  and  the 
impact  velocity,  v^.  From  the  usual  dimensional  considerations,  there¬ 
to  re,  the  solution  must  be  of  the  form 


where  r  and  z  are  cylindrical  coordinates  with  origin  at  the  center  of  im¬ 
pact  on  the  target  free  surface.  Thus,  the  variation  of  p^  causes  only  a 
proportional  change  in  the  densities  and  pressures  and  the  effect  of  vary¬ 
ing  Vq  is  only  to  change  the  time,  pressure,  and  velocity  scales  in  a 
manner  which  is  easily  determined  by  inspection  of  the  above  expressions. 
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Additional  simplification  of  the  solution  can  be  attained  if  we  introduce 
the  assumption  of  late -stage  equivalence.  The  assumption  is  suggested 
by  studies  of  solid-rolid  impact,  where  it  has  been  found  that  impacts  at 
different  velocities,  though  different  at  early  times,  rapidly  become  very 
nearly  the  same,  provided  the  relative  scale  of  the  two  impacts  is  such 
that  LqVq  is  constant.  The  analytical  statement  of  this  assumption  is  to 
require  that  the  solution  not  contain  L^,  v^  as  individual  parameters,  but 
only  in  the  combination  LqVq,  where  a  is  a  constant  to  be  determined.  It 
follows  directly  that  the  solution  must  approach  at  late  times  a  flow  express¬ 
ible  in  the  form 
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where 


T  =  tv0/L0,  R  =  r/L0,  Z  =  z/L0  . 

Thus  the  assumption  of  late -stage  equivalence  leads  one  to  a  similarity 
solution;  i.e.  ,  the  number  of  independent  variables  in  the  functions  is 
reduced  to  two  instead  of  the  original  three. 

The  reduction  to  two  variables  still  leaves  one,  however,  with  a  set 
of  partial  differential  equations,  and  the  flow  is  presumably  not  amenable 
to  analytical  solution  without  additional  approximations.  It  is  possible, 
however,  to  draw  certain  physical  conclusions  from  the  form  of  the  late- 
stage  solution; 
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1.  The  flow  configuration  is  asymptotically  self-similar.  The 
compression  p/p q,  for  example,  is  the  same  at  corresponding 
points  (equal  Z«+1/t“,  R«+1/t«).  Pressures  at  corresponding 

_  2  /  #y 

points  are  in  the  ratio  (Z^/Z^)  ,  where  Z^/Z^  is  the  ratio  of 

the  linear  dimensions  of  the  configuration  at  the  two  times. 
Similarly,  velocities  at  corresponding  points  drop  off  as  distance 
to  the  -(1/a)  power  as  the  disturbance  enlarges. 

2.  The  dimensions  of  the  disturbance  increase  as  ta!/(°!+l)j 

3.  The  mass  velocity  at  the  shock  front  and  the  shock  speed  attenuate 

with  penetration  distance  as  Z  or  t  The  shock  pres- 

..  *  7-2/a  .-2/(a+l) 

sure  attenuates  as  Z  or  t 

4.  The  total  positive  axial  and  total  positive  radial  momenta  increase 

as  t^  1)/(Q!+1)>  This  result  follows  by  noting  that  distances  in¬ 
crease  as  mass  engulfed  therefore  goes  as 

and  velocities  decrease  as  t  1/te+l)^  so  that  the  momentum  within 
corresponding  elements  of  volume  varies  as  t^  ^ 

The  Eulerian  code  has  been  used  to  treat  numerically  the  problem  of 

a  right  circular  cylinder  of  unit  aspect  ratio  impacting  a  semi-infinite 

target.  The  value  3/2  was  chosen  for  y  because  solids  at  extremely  high 

shock  pressures  (~  1,000  Mbars)  have  Hugoniot  curves  and  expansion 

isentropes  which  approximate  those  of  a  y  =  3/2  ideal  gas.  The  values 
6  3 

Vq  =  10  cm/ sec,  p^  =  1  g / cm  ,  and  a  projectile  diameter  of  3.  6  cm  were 
used  to  complete  the  problem  specification. 

A  point  of  particular  interest  concerning  the  computed  solution  is 
whether  or  not  the  flow  becomes  self-similar  as  indicated  by  the  assump¬ 
tion  of  late-stage  equivalence;  after  an  early  phase  (two  or  three  times 
the  time  required  for  the  initial  projectile  to  move  its  own  length),  the 
flow  geometry  tends  to  a  configuration  which  does  indeed  remain  self¬ 
similar  as  it  enlarges. 

In  this  late-stage  configuration,  the  shock  wave  closely  approximates 
a  spherical  shape,  with  the  (moving)  center  of  the  sphere  located  at  about 
one-third  of  the  shock  penetration  distance.  The  target  material  which 
has  been  engulfed  by  the  shock  is  largely  compressed  into  a  thin  shell 
near  the  shock  front.  Progressive  tages  of  this  flow  can  be  visualized 
from  the  velocity  and  mass  flux  plots  reproduced  as  Fig.  8. 


fields  for  the  axisymmetric  gas  problem  (arrows  indicate  vector 
lengths  for  T?  =  10"  cm/sec  and  (p"u/Pq)  =  10^  cm/sec) 


t  ■  9.75  ft  SEC 


t  ■  29  ft  SEC 

IX 10® 

t  *29  ft  SEC 


IX 10® 


Fig.  8  eont'd- -Velocity  (left)  and  mass  flux  or  momentum  (right) 
vector  fields  for  the  axisymmctric  gas  problem  (arrows  indicate 
ve'tor  lengths  for  tJ  =  10b  cm/sec  and  {pu/ Pq)  =  10°  cm/sec) 
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Tests  of  other  flow  properties  characteristic  of  the  similarity  solution, 
and  also  determinations  of  a,  are  possible  from  plots  of  the  integrated 
momenta  and  the  shock  attenuation.  Since  the  momenta  must  vary  as 
t(  3a  l)/(a+l)  ^  a  s  g  umpt  ion  of  late-stage  equivalence  is  valid,  one  can 

examine  the  log-log  plots  of  integrated  radial  and  axial  momenta  versus 
time  to  (a)  test  if  the  plots  become  linear  at  late  times,  (b)  test  whether 
the  late -time  slopes  of  the  two  momenta  plots  are  in  agreement,  and  (c) 
determine  a  from  the  slopes.  The  momentum  curves  are  plotted  as  Fig.  9, 
where  it  is  seen  that  the  plots  are  quite  linear  at  late  times  and  that  the 
slope  (3a  -  l)/(a  +  1)  is  in  good  agreement  with  both  plots  if  one  chooses 
a  =  0.  59.  (The  slope  for  a  =  1/  3  and  that  for  a  =  2/  3  are  indicated  on  the 
graph  for  comparison.  )  Figure  10  shows  log-log  plots  of  shock  penetra¬ 
tion  and  shock  strength  versus  time,  both  for  the  axis  of  symmetry.  The 
penetration  and  attenuation  slopes  are  in  reasonable  agreement  with  the 
values  a/(a  +  1)  and  -l/(or  +  1)  using  a.  =  0.  59.  (A  slightly  lower  a,  about 
0.  55,  would  correlate  the  results  somewhat  better;  but  the  value  of  a 
from  the  momentum  curves  is  considerably  more  accurate.) 

Finally,  velocity,  density,  and  pressure  profiles  at  various  times  can 
be  compared  to  test  whether  the  flow  remains  self- similar .  The  results 
of  a  representative  comparison- -the  velocity  distributions  along  45°  for 
times  different  by  a  factor  of  three  - -are  shown  as  Fig.  11.  It  is  seen  that 
if  the  curve  for  the  early  time  is  scaled  to  a  later  time  by  similarity  re¬ 
lations  Z~t  and  v~t  the  agreement  is  good.  The  value 

a  =  0.  59  was  used. 

The  results  of  all  comparisons  are,  then,  compatible  with  the  assump¬ 
tion  that  the  flow  has  tended  to  the  self-similar  form  indicated  by  late-stage 
equivalence.  It  is  also  of  interest  to  note  that  the  value  a  =  0.59  is  in  good 
agreement  with  the  values  obtained  for  condensed  materials,  despite  sub¬ 
stantial  equation-of- state  differences  and  formal  differences  in  the  method 
of  analysis  of  the  numerical  results. 

Using  the  assumption  that  the  flow  is  self-similiar,  Rae  and  Kirchner^ 
have  reported  some  very  interesting  analytical  progress  on  the  problem 
of  axisymmetric  ideal -gas  impact.  By  considering  the  flow  equations  as 
specialized  to  the  axis  of  symmetry  and  by  requiring  that  the  solution 
pass  smoothly  through  a  singularity  in  the  differential  equations,  they 


33 


*o  *o 

O3S/W0-9)  MfUNSWOW 


Fig.  9--Radial  and  axial  momenta  for  the  y  =  1.  5  ideal-gas  impact;  the  slopes  necessary  for 
late-stage  equivalence  on  an  Of  =  1/3  and  a  =  2/3  basis  are  indicated  at  the  upper  right 
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Velocity  versus  distance  alonp  line  45°  from  the  axis  of  symmetry;  the  scaled  curve  is 
determined  from  the  e?  rlier  data  and  the  similarity  relations 
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find  that  the  exponent  N  in  the  similarity  variable  r/t^  can  be  evaluated  to 
a  good  approximation.  Values  N  =  0.  37  to  0.  39  are  obtained,  as  com¬ 
pared  to  N  =  0.40  for  the  spherical-blast  problem. 

A  consequence  of  the  late  -  stage -equivalence  interpretation,  on  the 
other  hand,  is  a  relation  between  a  and  the  similarity  exponent  N.  Speci¬ 
fically,  comparison  with  the  present  similarity  variable  /rtt  shows 

that  N  =  a/(a  +  1).  Using  the  present  result,  a  =  0.  59,  gives  N  =  0.  37. 

Rae  gives  N  the  value  0.  375  for  y  =  3  and  the  value  0.  367  for  y  =  2.  Hence, 
the  value  for  y  =  3/2  is  found  by  extrapolation  to  be  0.  363. 

The  above  agreement  is,  of  course,  quite  gratifying.  The  numerical 

(9) 

integrations  and  the  analytical  treatment  by  Rae  and  Kirchner'  '  are  based 
on  entirely  different  approximations,  so  that  the  agreement  is  a  significant 
confirmation  of  both  results. 


The  principal  approximations,  in  addition  to  assuming  self-similarity, 
are:  The  shock  wave  in  the  target  is  a  hemisphere  and  an  estimate  is  made 
for  the  magnitude  of  d2p/d9Z  at  the  singular  point.  The  last  approximation, 
and  the  fact  that  N  is  found  to  be  a  weak  function  of  y,  accounts  for  the 
small  variation  in  N. 


7.  HYDRODYNAMICS  OF  SOLID-SOLID  IMPACT 


7.1.  DIMENSIONAL  ANALYSIS  AND  LATE -STAGE  EQUIVALENCE 

Previous  comparisons^^  showed  that  the  solid  equations  of  state 

appropriate  to  impact  work  can  be  correlated  by  assuming,  for  simplicity, 

that  two  properties,  the  initial  density  p  and  the  initial  sound  speed 

1/2  0 

c_  =  (dP /bp)  ,  suffice  to  characterize  a  given  material.  Specifically, 

u  s 

the  individual  equations  of  state  P  =  f(p,  E)  of  the  various  metals  are 

brought  into  reasonable  agreement  when  expressed  in  the  dimensionless 

2  2 
variables  p  =  P/pQcQ,  T?  =  p/ pQ,  and  e  =  E/cQ. 

For  fixed  initial  geometry,  say  that  of  a  sphere  or  cylinder  on  a 
thick  plate,  the  problem  parameters  for  a  like-material  impact  are  there¬ 
fore  the  material  properties  ,  c^,  the  impact  velocity,  v^,  and  a  charac¬ 
teristic  length,  such  as  the  projectile  diameter,  LQ .  From  dimensional 
considerations  the  solution  must  be  expressible  in  the  form 
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The  additional  assumption  that  two  impacts  be  late-stage  equivalent  on  the 
basis  of  equal  L^v®,  rather  than  depending  individually  on  the  parameters 
Lq,  Vq,  further  restricts  the  form  of  the  solution.  It  can  be  shown  that 
under  this  assumption,  the  solution  must  be  expressible  as 
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where  r  =  tUg/Lg,  R  =  r/L Q,  Z  =  z/Lq. 

From  the  first  form  of  the  solution,  one  can  conclude,  as  in  Ref.  1, 
that  the  solutions  for  two  geometrically  similar  impacts  of  the  same  vq/cq 
but  involving  different  materials  must  be  the  same  when  expressed  in 
dimensionless  variables. 

Additional  physical  conclusions,  similar  to  those  obtained  for  gases, 
do  not  appear  to  be  forthcoming.  It  is  not  possible,  for  example,  to  pre¬ 
dict  the  effects  of  impact  velocity  variation  from  the  dimensional  analysis 
results  because  of  the  (significant)  functional  dependence  on  v^/c^. 

Further,  the  late- stage  equivalent  form  of  the  solution,  though  simplified, 
still  involves  a  general  functional  dependence  on  three  (not  two)  independent 
variables.  Thus,  the  property  of  late- stage  equivalence,  although  verified 
by  the  numerical  solutions,  does  not  imply  a  similarity  solution  in  the  case 
of  solid-solid  impact. 

It.  is  apparently  necessary  to  treat  numerically  a  series  of  impact  prob¬ 
lems  in  order  to  determine  the  effects  of  such  variations  as  impact  velocity, 
the  shape  of  the  projectile,  or  the  equation  of  state  of  the  material.  To 
explore  these  and  similar  effects,  some  forty  impact  interactions  have 
been  computed  and  the  results  have  been  analyzed.  In  the  following,  the 
principal  conclusions  from  each  study  are  summarized  and  results  are 
presented  to  illustrate  specific  points. 
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7.2.  GENERAL  COMMENTS  ON  A  TYPICAL.  PROBLEM 

A  typical  interaction  in  the  present  series  is  that  of  a  10^  cm/sec 

cylindrical  lead  projectile  on  a  lead  target  plate,  the  projectile  being  3.  6 

cm  long  by  1.8  cm  radius.  The  continuous  Eulerian  code  was  used  and 

space  resolution  was  initially  defined  by  square  cells  0.  1  cm  on  a  side. 

Each  time  the  disturbance  reached  the  top  or  right  side  of  the  4300-cell 

grid,  the  cell  dimensions  were  doubled  by  means  of  a  rezoning  procedure 

and  additional  target  material  was  included.  The  problem  was  run  to 
2 

2.  1  5  x  10  |isec,  at  which  time  the  cell  dimensions  were  0. 8  cm,  the  axial 
shock  penetration  distance  was  approximately  56  cm,  and  the  shock  pres¬ 
sure  was  25  kbar. 

The  overall  features  of  the  flow  are  demonstrated  by  the  plots  of  the 
velocity  and  mass  flux  (or  momentum  per  unit  volume)  vector  fields  shown 
in  Figs.  12a  and  12b.  The  shock  front  in  the  metal  is  seen  to  rapidly  as¬ 
sume  an  approximately  hemispherical  shape.  The  motion  behind  the  shock 
wave  is  predominately  radial,  although  there  is  also  a  definite  tendency 
for  this  flow  to  diverge  away  from  the  axis  and  to  curl  outward  toward  the 
target  free  surface. 

Pressure  profiles  along  the  axis  of  symmetry  are  plotted  for  various 
times  in  Fig.  13.  Pressures,  especially  those  less  than  100  kbars  or  so 
in  a  metal,  are  the  least  precisely  defined  quantities  in  the  computation. 

In  particular,  oscillations  in  the  pressure  profile  are  apparent  behind  the 
shock  front  at  intermediate  times,  and  the  front  is  seen  to  be  quite 
smeared  out  at  late  times,  when  the  pressures  are  in  the  25-to  75-kbar 
range . 

The  profiles  give  a  reasonable  representation  of  the  overall  shock 
attenuation,  however,  and  are  also  in  very  good  agreement  with  the  ex¬ 
pected  plane-wave  shock  pressures  at  very  early  times.  In  general,  the 
Eulerian  results  on  cell-to-cell  variation  of  intensive  variables  have 
proved  to  be  more  reasonable  and  substantially  smoother  (see  Section  4) 
than  those  obtained  from  the  previous  particle -in-cell  calculations  and 
have  permitted  the  detailed  comparisons  of  solutions  from  different  prob¬ 
lems.  Other  plots  of  pressure  profiles  and  also  of  velocity  and  density 
profiles  are  presented  subsequently. 


Fig.  1  2b-- Velocity  (left-)  and  mass  flux  {right) 
vector  fields  for  the  10b  cm/ sec  lead  impact 
(arrows  indicate  vector  lengUi  for  Tf  =  10^ 
cm/  sec  and  pu/  =  10"*  cm/sec) 
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Integrated  or  averaged  extensive  quantities,  such  as  the  total  kinetic 
and  internal  energies  or  integrated  total  momenta,  are  much  more  smoothly- 
varying  than  the  intensive  cell  quantities.  Furthermore,  changes  in  mass 
resolution,  space  resolution,  minor  variations  in  the  hydrodynamic  code, 
and  variations  in  the  material  equation  of  state  have  had  very  small  effects 
on  these  quantities.  For  these  reasons  the  integrated  variables  of  the 
motion,  especially  of  the  total  momenta,  have  been  used  wherever  possible 
for  the  precise  comparison  of  flows. 

7.3.  VELOCITY  LATE -STAGE  EQUIVALENCE 

A  series  of  like-material  impacts  has  been  computed  for  iron, 
aluminum,  lead,  and  polyethylene.  In  each  (6-problem)  series  the  impact 

velocity  was  varied  from  1 .  24  c_  to  80  c_,  where  c_  is  the  material  "sound 

,1/2  U  U  U 

speed,"  c.  =  (3P/dp)  .  Thus,  the  velocities  considered  range  from 

those  attainable  experimentally  to  velocities  well  in  excess  of  ones  of 
interest  in  meteoroid  impact.  The  initial  geometry  for  each  interaction 
was  that  of  a  right  circular  cylinder  of  equal  length  and  diameter  inter¬ 
acting  with  a  thick  target  plate  of  the  same  material.  Parameters  defining 
the  various  problems  are  listed  in  Table  2. 


Table  2 

IMPACT  VELOCITIES  AND  PROBLEM  PARAMETERS  FOR  THE 
SERIES  ON  CYLINDRICAL  PROJECTILE  IMPACTS 


Iron 

Lead 

Aluminum 

Polyethylene 

Plastic 

vo/co 

Series 

(cm/  sec) 

(cm/ sec) 

(cm/ sec) 

(cm/  sec) 

A 

5  x  105 

2.  5  x  105 

6. 5  x  10 5 

3.  5  x  105 

1. 24 

B 

io6 

5  x  105 

1. 3  x  106 

7x  i°5 

2.  48 

A 

2  x  106 

106 

2. 6  x  106 

1. 4  x  106 

4.  96 

B 

4  x  106 

2  x  106 

5. 2  x  106 

2.  8  x  106 

9.  92 

A 

1. 6  x  107 

8  x  106 

2.  1  x  107 

1. 14  x  107 

40 

B 

3.  2  x  107 

1. 6  x  107 

4. 2  x  107 
_ 1 

2.  28  x  107 

80 

Note:  For  the  A  series,  projectile  length  and  diameter  were  3.6  cm 
and  cell  dimensions  prior  to  first  rezone  were  0.  1  cm;  for  the  B  series, 
the  projectile  dimensions  were  initially  2.4  cm  and  cell  size  was  the 
same  as  in  the  A  series.  Exceptions  are  the  first  four  iron  problems, 
which  were  somewhat  more  coarsely  zoned. 
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For  each  material  the  problems  were  planned  and  analyzed  in  pairs 
(1. 24  and  2. 48  c^,  4.  96  and  9. 92  c^,  40  and  80  c^).  The  purpose  was  to 
assure  comparable  space  resolution  at  late  times  so  that  comparison  of 
the  two  flows  for  late -stage  equivalence  could  be  made  without  possible 
bias  due  to  unequal  space  resolution.  In  each  pair,  the  twofold  increase 
in  impact  velocity  was  offset  by  reducing  the  projectile  dimensions  to  2/3 
the  former  values.  Thus,  the  problems  were  run  on  a  basis  intermediate 
between  that  given  by  keeping  projectile  kinetic  energy  or  projectile 
momentum  constant.  More  precisely,  the  indicated  change  in  velocities 
and  dimensions  corresponds  to  a  =  0.  585.  The  close  agreement  seen 
below  between  the  as-run  problems  (prior  to  changing  the  relative  size 
of  two  problems  by  scaling)  therefore  corresponds  to  values  of  a  which 
are  very  near  to  0.  585. 

For  the  mid-  and  high-velocity  pairs,  the  detailed  comparisons  of 
late  stages  of  the  interactions  show  very  good  agreement,  in  both  the 
plots  of  pressure,  density,  and  velocity  profiles  from  individual  cell 
quantities  and  in  the  integrated  momentum  curves.  Typical  results  of 
these  comparisons  are  reproduced  as  Figs.  14,  15,  and  16.  The  very 
good  late-stage  agreement  of  the  flows  is  a  gratifying  confirmaion  of  the 
equivalence  first  reported  from  the  particle -in-cell  calculations.^  The 
present  results  considerably  strengthen  the  earlier  basis  for  claiming 
equivalence  because  the  smoother  solutions  make  possible  quantitative 
comparisons  of  the  intensive  cell  variables,  whereas  the  previous  particle- 
in-cell  results  necessarily  relied  primarily  on  integrated  quantities. 

For  the  lowest  velocities  considered  (i.e.  ,  5  x  10^  to  10^  cm/sec  in 
iron,  2.  5  x  10  to  5  x  10  cm/fec  for  lead,  etc.),  late-stage  equivalence 
is  less  exact.  The  radial  and  axial  momentum  curves  are  again  brought 
into  good  agreement  by  a  single  value  of  ol,  but  there  persists  to  late 
times  some  tendency  for  the  higher  velocity  impact  to  be  associated  with 
a  more  intense  but  narrower  pulse.  Thus,  in  the  iron  pair,  the  shock  for 
the  10^  cm/ sec  impact  has  attenuated  to  250  kbar  while  that  for  the  low- 
velocity  impact  is  down  to  200  kbar  at  the  same  radius;  for  lead,  an 
analogous  comparison  is  90  and  70  kbar.  It  seems  apparent  that  the 
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Fig.  14.  Convergence  to  late-stage  equivalence 
of  two  impact  flows.  Pressure  versus  distance 
at  three  times  and  o/pq  versus  distance  at  the 
latest  time.  Plots  are  from  the  8  x  10°  cm/sec 
and  1.6  x  10^  cm/sec  lead  impacts.  Distance  is 
measured  along  a  line  45°  from  the  axis  of  sym¬ 
metry.  with  origin  at  the  center  of  impact  on  the 
target  surface 
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iparison  of  axial  pressure  profiles  for  the  2  x  10°  cm/sec  and  4  x  10°  cm/sec 
iron  impacts,  showing  late-stage  agreement 
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failure  to  achieve  precise  equivalence  is  characteristic  of  low  velocities 
and  would  become  more  pronounced  if  the  velocity  were  further  reduced. 

Values  of  Ct  were  obtained  from  the  various  comparisons  by  scaling 
the  relative  size  of  the  problems  until  the  best  agreement  is  obtained 
between  the  two  solutions.  (The  most  precise  scale  factors  are  those 
given  by  requiring  that  the  axial  and  radial  momenta  agree.  )  is  then 
calculated  from  L^v^*  =  L^v®,  where  L^  and  are  the  characteristic 
linear  dimensions  after  scaling  and  v^,  v^  are  the  impact  velocities. 

The  results  show  a  remarkable  uniformity  in  ex  for  all  twelve  comparisons. 
Individual  values  for  the  eight  higher-velocity  pairs  lie  between  0.  56  and 
0.  59  and  the  low-velocity  comparisons  indicate  values  of  ex  very  near  0.60. 
Precisions  on  the  individual  determinations  are  typically  3%  to  4%. 

The  present  results  are  in  substantial  agreement  with  previous  values 
(ex  0.62)  which  were  obtained  for  iron  from  particle-in-cell  calculations. 
The  offset  of  approximately  0.  0^x  is  due  to  a  more  careful  interpretation  of 
the  momentum  curves,  and  the  present  value  «  ='  0.  58  should  be  considered 
more  accurate. 

7.  4.  MATERIAL  EQUIVALENCE 

As  was  noted  previously,  the  c^  equaHon-of- state  correlation 
implies  that  geometrically  similar  impact  problems  of  equal  dimensionless 
impact  velocity,  Vq/c q>  should  have  the  same  solution  in  dimensionless 
variables.  This  property  of  impact  hydrodynamics  is  of  some  interest 
because  it  enables  one  to  generalize  computed  (or  experimental)  hydro- 
dynamic  results  learned  from  work  on  only  one  or  two  representative 
materials . 

Comparisons  to  test  the  proposed  equivalence  for  iron,  lead,  aluminum 
and  plastic  interactions  are  possible  from  any  pair  of  problems  listed  in  a 
given  line  of  Table  2.  Dimensionless  plots  of  the  momentum  curves  and 
representative  pressure  profiles  have  been  made  for  each  of  the  possible 
comparisons  and  typical  results  are  reproduced  as  Figs.  17  through  20. 

The  general  agreement  is  seen  to  be  very  good.  It  may  be  noted  that  the 
results  for  polyethylene  correlate  somewhat  less  satisfactorily  than  those 
for  the  metals,  as  should  be  expected  from  the  fact  that  polyethylene  was  a 
nonconforming  member  of  the  underlying  equation- of- state  correlation. 
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Fig- 


19- -Comparison  of  momenta  versus  time  in  dimensionless  variables 
for  polyethylene,  aluminum,  and  lead  impacts  of  equal  Vq/cq 
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Additional  tests  of  material  equivalence  scaling,  based  on  results 
from  particle-in-cell  calculations,  are  described  in  Ref.  1.  The  com¬ 
parisons  given  there  are  limited  to  integrated  or  averaged  quantities  but 
are  in  good  accord  with  the  conclusions  given  here.  The  above-reported 
agreement  in  the  values  of  a  for  the  various  materials  is  also  to  be  expected 
from  the  equation-of- state  correlation. 

7.  5.  PROJECTILE -SHAPE  EFFECTS 

A  series  of  twelve  impact  interactions  were  studied  in  order  to  explore 
effects  of  projectile  shape.  Four  shapes  were  considered  for  each  of  three 
impact  velocities: 

1.  Right  circular  cylinder,  equal  length  and  diameter, 

2.  Sphere, 

3.  Rod,  three  diameters  long, 

4.  Disc,  three  thicknesses  in  width. 

The  velocities  were  1. 24  Cq,  4.  96  cQ,  and  40  cQ  (corresponding  to  three 
of  the  cylindrical  projectile  impacts  listed  in  Table  2)  and  the  projectile 
masses  were  the  same  in  all  problems.  The  impacting  materials  were 
taken  to  be  lead,  although  the  results  apply  also  to  other  solids  by  virtue 
of  the  material  equivalence  studies. 

In  the  analysis  of  the  computed  results,  particular  attention  was  given 
to  determining  whether  the  different  projectile  shapes  lead  to  a  late-stage 
equivalence  and  to  evaluating  the  relative  effectiveness  ox  the  various  shapes. 

The  simplest  and  most  definitive  results  were  obtained  in  the  com¬ 
parisons  of  spherical  and  cylindrical  (shape  1, above)  projectiles.  For 
each  of  the  three  impact  velocities,  these  projectiles  lead  to  practically 
identical  late-stage  flows,  both  in  the  total  axial  and  radial  momenta  and 
in  the  comparison  of  instantaneous  pressure,  density,  and  velocity  plots. 
Substantial  differences  are  limited  to  the  very  early  times  in  which  the 
projectile  is  embedding  itself  in  the  target  plate.  Finally,  the  relative 
scale  of  the  two  flows  is  such  that  the  projectiles  are  equivalent  on  an 
equal-mass  basis. 

The  equivalence  of  spheres  and  cylinders  has  been  of  some  concern 
in  previous  investigations,  with  opinion  being  divided  between  equivalence 
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on  the  basis  of  equal  projectile  mass  and  the  substantially  different  con¬ 
clusion  that  a  sphere  of  given  diameter  is  equivalent  to  a  unit  aspect  ratio 
cylinder  of  that  same  diameter.  Theoretical  or  computational  work  on  the 
problem  has  not  been  previously  reported,  however,  and  the  present  result 
is  apparently  in  accord  with  more  recent  experimental  findings. 

Late-stage  equivalence  is  somewhat  less  exact  for  the  disc  and  rod 
projectile  geometries.  At  the  two  higher  velocities,  however,  the  flows 
are  in  reasonably  good  general  agreement  after  a  shock  penetration 
distance  which  is  large  compared  to  initial  projectile  dimensions.  This 
is  demonstrated  by  the  vector  field  plots  for  the  4.  96  c^  (10^  cm/sec  for 
lead)  impact,  as  seen  in  Figs.  2,1  and  22.  A  close  comparison  of  these 
results,  however,  shows  that  even  at  late  times  the  disturbance  for  the 
rod  impact  is  slightly  more  elongated  in  the  axial  direction  than  that  for 
either  the  disc  or  the  sphere.  Except  for  this  characteristic  of  the  rod 
impacts,  the  late-stage  agreement  found  in  the  high-velocity  comparisons 
(10^  and  8  x  10^  cm/ sec)  is  quite  good.  Precise  comparisons  of  the 
momentum  curves  and  the  associated  scale  factors  show  that  the  rod 
impact  is  equivalent  to  that  for  the  sphere  or  cylinder  on  an  projectile 
equal  mass  basis.  The  disc  is  slightly,  but  perceptibly,  less  effective. 

At  both  velocities,  it  would  require  an  increase  of  some  15%  in  disc  mass 
for  this  projectile  to  be  as  effective  as  the  others. 

The  above  rod  and  disc  projectiles  were  purposely  chosen  to  have 
moderate  aspect  ratios,  with  the  objective  of  testing  the  limits  of  a  late- 
stage  equivalence  to  show  independence  of  projectile  shape.  Substantially 
more  extreme  aspect  ratios  will  certainly  lead  to  flows  which  remain 
qualitatively  different  throughout  the  hydrodynamic  regime.  Such  dissimilar 
flows  are,  in  fact,  also  evident  for  the  present  aspect  ratios  at  the  lowest 

5 

velocity  (2.  5  x  10  )  studied.  The  target-plate  disturbance  due  to  the  rod, 
for  example,  was  much  more  concentrated  near  the  axis  than  that  for  the 
disc,  and  the  impact  configurations  are  expected  to  remain  dissimilar 
indefinitely. 

Analytical  results  can  be  cited  for  the  limiting  extremes  of  thin-disc 
and  long-rod  geometries.  For  disc  projectiles  of  sufficiently  extreme 
aspect  ratio,  impact  disturbances  will  become  equivalent  as  essentially 
plane  waves.  The  slab  impact  considerations  for  an  ideal  gas  of  y  =  1.5, 
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Fig.  22--Comparison  of  the  mass  flux  on  momentum  distributions 
for  the  sphere,  rod,  and  disc  impacts 

in  particular,  indicates  equivalence  on  a  basis  which  is  approximately 
midway  between  that  for  equal-energy  and  equal-momentum  projectiles. 
Because  of  equation-of-state  similarities,  the  result  should  apply  also  to 
solid  impact  for  sufficiently  high  velocities  and  sufficiently  extreme  disc 
aspect  ratio.  It  may  be  that  the  relative  ineffectiveness  of  the  aspect 
ratio  three  disc  problems  is  indicative  of  the  fact  that  a  should  decrease 
to  about  0.  51  in  this  limit. 

The  classical  jet-penetration  theory  applicable  to  long  rod  geometries 
is  sufficiently  well  known  that  review  here  is  uncessary.  It  may  be  noted, 
however,  that  recent  applications^  of  a  jet-penetration  model,  by  Allison 
and  Vitali  of  the  Ballistic  Research  Laboratories,  have  led  to  an  extremely 
encouraging  correlation  with  experimental  shock-wave  observations  of  a 

5 

9  x  10  cm/sec  aluminum  jet  on  a  plexiglass  target.  Such  considerations 
are,  indeed,  quite  applicable  to  the  early  stages  of  the  rod  impacts  dis¬ 
cussed  above . 

7.6.  DENSITY-EFFECT  STUDIES 

In  most  applications,  the  projectile  and  target  are  different  materials, 
and  it  is  desirable  to  extend  present  considerations  to  include  such  cases. 
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Dimensional  considerations  suggest  that  one  should  introduce  two  new 
parameters  to  characterize  unlike-material  impacts--namely,  the  density 
and  sound  speed  of  the  projectile  in  addition  to  those  of  the  target.  It 
seems  reasonable  to  assume,  however,  that  variations  in  projectile 
density  are  more  significant  than  those  of  sound  speed,  and  the  subject  of 
unlike-material  impact  is  usually  equated  to  that  of  density-effect  studies. 

Several  comparisons  have  been  made  to  evaluate  possible  effects  due 
to  changing  projectile  density.  These  include: 

1.  Three  ideal-gas  impacts,  y  =  1.5,  were  computed  in  which  the 

3 

initial  projectile  density  was  (1.5)  =  3.  375  times  the  target  plate  density, 

.  3 

equal  to  the  target  density,  and  (1.5)  times  the  target  density;  the  initial 
masses  of  the  three  projectiles  were  the  same.  The  problems  were  run 
using  the  continuous  Eulerian  code.  The  calculations  were  carried  to  a 
time  about  125-fold  greater  than  the  time  required  for  the  pellet  to  bury 
itself  in  the  target,  and  detailed  cross  comparisons  of  the  flows  were  made 
in  order  to  evaluate  late-stage  equivalence  and  the  relative  effectiveness  of 
the  three  densities.  When  consideration  is  given  to  the  fact  that  the  total 
variation  in  projectile  density  is  by  a  factor  of  11.4,  the  flows  are  in 
remarkably  good  agreement.  Shock  penetration  distances  at  late  times 
agree  to  10%,  with  the  maximum  discrepancy  occuring  on  the  axis  of 
symmetry.  Axial  momenta  in  the  three  cases  agree  within  5%,  and  the 
spread  in  radial  momentum  is  20%  at  late  times.  The  direction  of  indicated 
discrepancies  is  such  that  the  high-density  projectile  is  a  more  effective 
projectile,  but  the  effect  is  a  negligibly  small  one. 

2.  Two  calculations  were  made  using  special  equations  of  state  in 

order  to  investigate  possible  projectile  density  effects  for  metal-metal 

impacts  at  velocities  around  10^  cm/sec.  The  target  material  was  normal 

lead  and  the  projectiles  were  taken  to  have  a  modified  lead  equation  of 

3  -  3 

state  for  which  the  densities  were  respectively  (1.  5)  and  (1.  5)  times 
those  of  normal  lead.  This  procedure  was  used  to  provide  a  comparison  in 
which  the  only  equation- of- state  variation  was  the  density  ratio  for  the 
target  and  projectile.  The  calculations  were  carried  out  using  the  particle- 
in-cell  code,  and  a  relatively  high  mass  resolution  of  36  particles  per  cell 
was  used  in  order  to  facilitate  the  comparison  of  cell  quantities.  The  two 
problems  were  run  with  equal  initial  projectile  masses.  The  remarks 
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made  above  for  gas  impacts  apply  almost  verbatim  to  the  present  com¬ 
parison.  Axial  momenta  agree  well  and  the  total  spread  in  radial 
momentum  at  late  times  is  only  24%.  The  effects  of  the  11.4-fold 
variation  in  projectile  density,  with  projectile  mass  kept  constant,  are 
quite  small;  within  20%  in  projectile  mass,  the  two  projectiles  are 
equivalent  on  an  equal-mass  basis. 

3.  Three  impacts  were  computed  in  which  tungsten,  iron,  and 
beryllium  projectiles  interacted  with  iron  targets  at  a  velocity  of  4  x  10^ 
cm/sec.  The  particle -in-cell  code  was  used;  the  interactions  were  again 
found  to  be  approximately  equivalent  at  late  times  for  equal  initial  pro¬ 
jectile  masses. 

Substantial  premiums  on  projectile  density  have  occasionally  been 
reported  from  empirical  correlations^^’  of  experimental  impact  data. 
The  problem  has  not  been  the  subject  of  previous  theoretical  treatment, 
however,  and  the  present  results  are  believed  to  be  sufficient  to  show 
that  such  effects  do  not  persist  to  velocities  of  10^*  cm/sec  and  greater. 


8.  OVERALL,  EFFECTS  OF  SOLID-SOLID  IMPACT 


8.  1.  PENETRATION  RELATION 

The  close  agreement  of  two  impact  flows  during  the  hydrodynamic 
phase  of  the  interaction  is  useful  when  one  observes  that  the  subsequent 
strength-dependent  phase  of  the  interaction  will  also  be  the  same.  It 
follows  that  the  craters  in  the  two  cases  must  agree  in  shape  and  in  size. 
This  result  and  the  equivalence  observations  reported  in  Section  7  can  be 
summarized  by  the  penetration  relation 
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Here  p  is  the  crater  depth  or  some  other  characteristic  crater  dimension, 
d  is  the  diameter  of  the  projectile  if  the  mass  were  considered  to  be  a 
sphere,  and  ct  are  target  density  and  sound  speed,  p  and  v  are  pro¬ 
jectile  density  and  velocity,  and  K  is  a  constant  to  be  determined  from  a 
single  impact  experiment.  The  quantity  K  is  a  constant  of  the  target 
material,  depending  on  its  strength. 

From  the  results  given  in  Section  7,  one  can  expect  the  penetration 
relation  to  be  reliable  for  solid  impacts  in  which  (1)  the  impact  velocity 
is  greater  than  1—2  times  target  sound  speed;  (2)  the  projectile  aspect 
ratio  is  within  a  factor  of  3  of  unity;  and  (3)  the  projectile  density  is 
within  a  factor  of  about  3  of  the  target  density,  which  may  be  too  con¬ 
servative,  but  additional  calculations  would  be  required  to  confirm  this. 


8.2.  MOMENTUM  MULTIPLICATION 

Another  interesting  overall  effect  of  a  hypervelocity  impact  is  that  the 
momentum  transmitted  to  the  target  may  be  many  times  the  momentum  of 
the  incident  projectile.  The  dependence  of  this  momentum  upon  impact 
velocity  is  easily  determined  by  recalling  from  Section  7  that  the  late- 
stage  momentum  depends  only  on  MqVq0^  rather  than  on  M^,  v^  individually. 
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That  is, 

I  =  f(M0v03“). 

On  the  other  hand,  simple  linear  scaling  shows  that  impulse  must  be 
directly  proportional  to  projectile  mass  Mq,  so  that  this  functional 
dependence  must  be 
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where  k  is  a  constant.  Finally,  momentum  multiplication,  defined  as  the 
ratio  of  total  target  impulse  to  initial  projectile  momentum,  satisfies 
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Hence,  momentum  multiplication  increases  as  the  0.  74  power  of  impact 
velocity,  when  the  value  O'  -  0.  58  is  taken  into  account. 


8.  3.  TARGET -HARDNESS  EFFECTS 

"'here  has  been  some  interest  in  the  past  as  to  whether  target  hardness 
(known  to  play  an  important  role  in  limiting  crater  size  at  low  impact  veloci¬ 
ties)  persists  in  dominating  crater  size  as  the  velocity  is  increased  The 
affirmative  answer  to  this  question  leads  to  the  conclusion  that  a  reasonable 
estimate  of  crater  size  cannot  be  made  using  theoretical  hydrodynamic 
results  only. 

An  important  experimental  study  of  this  question  was  made  recently  by 

(13)  5 

J.  H.  Kineke  and  L.  G.  Richards.  *  '  Velocities  of  9.  7  x  10  cm/sec  and 
1.  55  x  10^  cm/ sec  were  attained  for  aluminum  and  beryllium  projectiles, 
respectively,  both  of  which  were  fired  into  soft  and  hard  aluminum  targets. 
The  ratio  of  crater  volumes  in  the  two  targets  (same  impacting  projectile) 
was  found  to  be  a  constant,  about  3,  for  the  entire  range  of  impact  veloci¬ 
ties  from  a  few  times  10^  cm/sec  to  the  maximum  values  given  above. 

Thus,  target- strength  effects  continue  to  influence  crater  dimensions  up 
to  the  highest  impact  velocities  attained. 

The  Kineke -Richards  result  is  precisely  the  one  which  should  be 
expected  from  late-stage  equivalence .  Consider,  for  example,  a  pair  of 
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impacts  (hard  and  soft  targets)  at  5  x  105  cm/sec.  The  pair  will  obviously 
be  identical  in  the  hydrodynamic  phase,  because  the  hydrodynamic  equa¬ 
tions  of  state  are  the  same.  This  pair  is  also  late-stage  equivalent  to 
another  pair  (hard  and  soft  targets)  at  1.  5  x  10^  cm/sec.  Hence,  all  four 
impacts  lead  to  the  same  late-stage  deformation  problem.  A  ratio  of  3  in 
crater  volumes  for  the  low-velocity  pair,  in  particular,  must  also  apply 
to  the  1.  5  x  10^  cm/sec  experiments.  Target- strength  effects  do  not 
become  less  important  as  impact  velocity  is  increased. 


9.  CONCLUDING  REMARKS 


9.1.  PRINCIPAL  CONCLUSIONS 

The  principal  conclusions  from  the  present  work  are  that 

1.  Crater  dimensions  increase  as  the  0.  58  power  of  the  impact 
velocity. 

2.  Changing  projectile  density  at  constant  projectile  mass  leaves 
unchanged  the  overall  effects  on  the  target. 

3.  Spherical  and  unit-aspect- ratio  cylindrical  projectiles  of  equal 
mass  have  equivalent  impact  effects  except  at  very  early  times. 
Projectile  aspect  ratios  as  extreme  as  3  are  unimportant  for 
impacts  at  five  times  target  sound  speed  and  greater. 

4.  Target  strength  plays  an  important  role  in  limiting  crater  size 
even  for  arbitrarily  high  velocities.  The  ratio  of  crater  dimen¬ 
sions  for  hard  and  soft  targets  is  independent  of  the  impact 
velocity. 

5.  The  ratio  of  total  target  impulse  to  initial  projectile  momentum 
increases  as  the  0.  74  power  of  the  impact  velocity. 

Conclusions  (1),  (4),  and  (5)  are  applicable  to  impact  velocities  which  are 
greater  than  about  twice  target  sound  speed,  and  (2)  pertains  to  impact 
velocities  greater  than  about  five  times  target  sound  speed. 

Conclusion  (1)  is  in  good  general  agreement  with  the  0.  62  power  de¬ 
pendence  reported  previously,^  The  basis  for  this  result  has,  however, 
been  considerably  strengthened  in  the  present  investigation:  The  hydro- 
dynamic  calculations  have  been  extended  to  much  later  times  and  to  pres¬ 
sures  (0.  1  Mbar  and  less)  which  are  typically  an  order  of  magnitude  lower. 
More  important,  the  smoothness  of  the  continuous  Eulerian  solutions  per¬ 
mits  one  to  make  detailed  comparisons  of  plots  of  the  intensive  variables 
made  from  listings  of  individual  cell  quantities.  Such  comparisons  (e.g.  , 
Fig.  14)  confirm  the  earlier  claim  that  flows  at  different  impact  velocities 
are  late-stage  equivalent. 


61 


62 


9.2.  EXPERIMENTAL -THEORETICAL  CORRELATIONS 

Results  of  the  present  investigation  appear  to  be  in  good  agreement 
with  the  main  body  of  experimental  data.  The  results  of  comparisons  with 
recent  experimental  findings  are  briefly  summarized  below: 

1.  Penetration  Exponents.  The  value  2/3  for  a  has  frequently  been 

quoted  from  summaries^  of  available  experimental  knowledge.  Energy 

(13) 

equivalence  (a  =  2/3)  was  also  reported  by  Kineke  and  Richards  '  '  from 

aluminum  data  extending  to  9.7  x  10  cm/sec.  Also  Halperson's  recent 

(15)  5 

1100F  data  for  aluminum,  '  extending  to  7  x  10  cm/s ec,  appear  to  us 
to  be  quite  reasonably  fitted  by  a  constant  exponent  a  =  1/2. 

The  present  theoretical  result  for  aluminum  is  that  a  =  0.58  for  veloci¬ 
ties  above  about  10^  cm/sec.  The  agreement  with  experiment  seems  very 
satisfactory,  especially  when  account  is  taken  of  the  fact  that  the  experi¬ 
mental  values  must  be  obtained  by  differentiating  penetration  data  and  ac¬ 
cordingly  apply  for  velocities  somewhat  less  than  the  upper  limit  of  the 
data. 

2.  Momentum  Multiplication.  Total  target  impulses  that  are  greater 
than  the  initial  pellet  momentum  have  been  observed  experimentally.  In 

a  previous  progress  report  (see  Ref.  20a)  it  was  shown  that  recent  data^**^ 
for  lead  were  compatible  with  conclusion  (5), above. 

3.  Target  Strength.  The  experimental  findings  by  Kineke  and  Richards 
were  discussed  in  Section  8  and  seen  to  be  in  agreement  with  the  theoretical 
result  on  persistence  of  strength  effects  at  high  velocities. 

9.3.  COMPARISON  WITH  OTHER  THEORETICAL  RESULTS 

The  first  application  of  hydrodynamic  codes  to  the  impact  problem  was 

( 17) 

the  well-known  work  by  Bjork.  Aluminum -aluminum  and  iron-iron  im- 

5  6 

pacts  were  treated  at  velocities  of  5.  5  x  10  cm/sec,  2x10  cm/ sec,  and 
7.2  x  10^  cm/sec. 

While  present  results  and  those  by  Bjork  are  apparently  in  reasonable 
agreement  for  the  so-called  hydrodynamic  phase,  an  important  difference 
in  the  two  approaches  arises  in  the  methods  used  to  determine  the  penetra¬ 
tion  relation.  Bjork  obtained  penetrations  by  estimating  the  stage  at  which 
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strength  effects  would  arrest  the  hydrodynamics.  From  three  such  esti¬ 
mates  for  each  metal  he  concluded  that  penetration  increases  as  the  1/3 
power  of  the  impact  velocity. 

/  io\ 

In  his  recent  discussion'  '  for  the  Sixth  Symposium  on  Hypervelocity 
Impact,  however,  Bjork  revised  his  earlier  penetrations  for  impact  veloci¬ 
ties  greater  than  about  2  x  10  cm/sec.  The  new  penetration  for  velocity 
7.  2  x  10^  cm/ sec  was  estimated  to  be  some  40%  higher  than  the  old  one. 
The  correction,  made  to  take  better  account  of  melting,  was  given  as  a 
rough  estimate.  This  estimated  correction  is  quite  important,  accounting 
for  about  one-half  the  discrepancy  between  Bjork's  old  data  and  our  pres¬ 
ent  result. 

It  seems  desirable  to  await  a  more  complete  evaluation  of  the  melting 
effect  before  attempting  to  decide  whether  a  significant  difference  exists 
in  the  two  approaches.  It  is  of  interest  to  note,  however,  that  the  remain¬ 
ing  discrepancy  is  substantially  less  than  that  for  the  extremes  of  energy 
equivalence,  a  =  2/3,  and  momentum  equivalence,  a  =  1/3.  These  Qf's, 
if  used  to  extrapolate  the  experimental  data  by  a  factor  of  10,  lead  to  pene¬ 
trations  (suitably  chosen  scale)  of  214  and  100  fora  =  2/3  and  a  =1/3, 
respectively.  On  this  same  scale  the  penetration  for  a  =  0.  58  is  174  and 
Bjork's  revised  penetration  is  140. 

It  should  be  remarked,  in  connection  with  melting,  that  the  method 
used  in  the  present  effort  does  not  require  a  knowledge  of  the  effects  of 
melting  or  of  shock  heating  on  target  strength.  The  approach  depends  on 
showing  merely  that  two  impacts  are  hydrodynamically  equivalent  prior 
to  the  onset  of  the  strength-dependent  deformation  phase. 

(19) 

Riney  has  very  recently  published  some  interesting  new  results. 

A  particle -in-cell  code  was  applied  to  aluminum-aluminum  and  lead-lead 
impacts,  and  the  conclusion  was  reached  that  impacts  are  late-stage 
equivalent  on  an  equal-energy,  or  a  -  2/3,  basis.  It  was  further  concluded 
that  target-strength  effects  do  not  vanish  as  impact  velocity  increases. 

The  general  agreement  between  Riney's  results  and  the  present  in¬ 
vestigation  is  very  encouraging.  Riney  does  not  claim  a  theoretical  basis 
for  exact  energy  equivalence,  however,  and  none  has  ever  been  given. 
Indeed,  our  most  important  general  conclusion  from  the  one-dimensional 
impact,  the  axisymmetric  ideal-gas  impact,  and  the  work  on  solid-solid 
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impact  is  that  impacts  are  equivalent  on  a  basis  intermediate  between  equal 
projectile  momentum  and  equal  projectile  energy.  For  the  important  case 
of  solid-solid  impact  at  velocities  in  the  meteoric  range,  energy  equivalence 
is  more  nearly  correct. 
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